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Abstract 



We calculate the tunneling conductance spectra of a ferromagnetic metal 
/ insulator / triplet superconductor from the reflection amplitudes using the 
Blonder-Tinkham-Klapwijk (BTK) formula. For the triplet superconductor, 
we assume one special p-wave order parameter, having line nodes, and two 
two dimensional /-wave order parameters with line nodes, breaking the time 
reversal symmetry. Also we examine nodeless pairing potentials. The evolu- 
tion of the spectra with the exchange potential depends solely on the topology 
of the gap. The weak Andreev reflection within the ferromagnet results in the 
suppression of the tunneling conductance and eliminates the resonances due 
to the anisotropy of the pairing potential. The tunneling spectra splits asym- 
metrically with respect to -E = under the influence of an external magnetic 
field. The results can be used to distinguish between the possible candidate 
pairing states of the superconductor Sr2Ru04. 
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I. INTRODUCTION 



The recent discovery of superconductivity in Sr2Ru04 has attracted much theoretical 
and experimental interest [|l| . The time reversal symmetry is broken for the superconductor 
Sr2Ru04, and the magnetic field is spontaneously induced as shown by /iSR experiment. 
0. The Knight-shift shows no change when passing through the superconducting state and 
is a clear indication for spin triplet pairing state, with a d- vector aligned to the z-a.xis. [Q. 
In addition the band structure calculations Q and de Haas van Alphen measurements 
show little dispersion along which is consistent with a two dimensional basis function on 
a cylindrical Fermi surface. Furthermore, the presence of a large residual density of states 
of quasiparticles inside the superconducting gap is evident from the linear temperature 
dependence of the nuclear spin lattice relaxation rate 1/Ti of ^°^Ru bellow OAK. Also 
specific heat measurements support the scenario of line nodes within the gap as in the high 
Tc cuprate superconductors 0. 

In the tunneling experiments between normal metals and superconductors Andreev re- 
flection process take place In the Andreev reflection process an electron incident, in 
the barrier with an energy bellow the gap can not drain off into the superconductor. It 
is instead reflected as a hole and a Cooper pair is transfered into the superconductor. In 
anisotropic high superconductors due to the sign change of the pair potential that the 
transmitted quasiparticles feel, zero energy states are formed, which are detected as peaks 
in the conductance spectra at -E = 0. Also in the presence of an imaginary s-wave 
component, which breaks the time reversal symmetry, the zero energy peak is shifted to the 
amplitude of the subdominant component. |[lT|jr^ 



The properties of the Andreev reflection are modified in the presence of an exchange 
field as in ferromagnet / insulator / superconductor junctions, since the retro-reflection of 
the Andreev reflection is broken in the ferromagnet. This phenomenon has been clarified 
both in s-wave and d-wave junctions and interesting aspects of the Andreev reflection have 
been revealed. [|T^-|T6[ Also the properties of ferromagnet / insulator / triplet superconductor 
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junctions have been studied where two types of pairing potentials are assumed for the triplet 



superconductor, i.e. the unitary and the non unitary with symmetry. [|T^ In the unitary 
case the conductance within the gap is reduced with the exchange interaction while in the 
non unitary it is not much influenced since the latter pairing state conserves spin. 

In order to identify the pairing state of Sr2Ru04 the Bogoliubov-de Gennes (BdG) equa- 
tions have been used to calculate the quasiparticle bound state wave function around non- 
magnetic impurities in unconventional superconductors. |]18| The characteristic patterns were 
distinguished for two proposed order parameters, (i) Eu, and (ii) Big x Eu- 

In this paper we will use the BdG equations to calculate the tunneling conductance of 
ferromagnet / triplet superconductor contacts, with a barrier of arbitrary strength between 
them, in terms of the probability amplitudes of Andreev and normal reflection. For the 
triplet superconductor we shall assume three possible pairing states of two dimensional 
order parameters, having line nodes within the Ru02 plane, which break the time reversal 



symmetry. The first two are the 2D /-wave states proposed by Hasegawa et al, having 
Big X Eu and B2g x E^ symmetry respectively. The other one is called nodal p-wave state 
and has been proposed by Dahm et al where the pairing potential has the form d{k) = 
Aoi;(sin(A;^a) +ism{kya)), with k^a = ttcosO, kya = vrsin^. This pairing symmetry has 
nodes as in the B2g x Eu case. Also we will consider two nodeless pairing states. One 
is the isotropic p-wave state and the other is the nodeless p-wave state initially proposed 
by K. Miyake and O. Narikiyo both breaking the time reversal symmetry. Generally 
the tunneling conductance is suppressed with the increase of the exchange interaction and 
the peaks are removed. This is due to the suppression of the Andreev reflection in the 
ferromagnet. For the nodal pairing states the linear dependence of the tunneling conductance 
with E is not much influenced. For the nodeless cases, the normalized conductance develops 
a constant value within the gap, which is suppressed as the exchange field gets larger. 
When the ferromagnet is a normal metal, the magnetic field splits the tunneling spectrum 
symmetrically around E = 0. The exchange field eliminates the negative branch of the 
tunneling spectra in the half metallic ferrogmagnetic limit. 



II. THEORY OF TUNNELING EFFECT 



For spin-triplet superconductors the wave functions describing the quasiparticles ^(r) 
are four-spinors in Nambu (particle-hole ® spin) space. Their particle and hole components 



are determined by the solutions of the BdG equations p3|,pl 



where, 



E^(r) = / dr'H(r,r')'^(r') 



(1) 



(2) 



H(r,r') 



He A 
A* -H, 



(3) 



A is the 2x2 triplet pairing matrix with elements of the form Ass{r, r'), and the spin index 
s=t, or s =|. He = He{r')6{r — r')ao, where do is the 2x2 unit matrix, and Heir) is 
the single-particle Hamiltonian which is given by Tieir) = —h^ Vr /2^e + y{f) ~ Ep, E is 
the energy measured from the Fermi energy Ep. For the pairing states we examine in Sec. 
Ill, the spin-up and spin-down components decouple. We will consider only triplet pairing 
states where A||(r,r') = A||(r,r') = 0, while A^{r,r') = Aj^|(r,r'). In that case the 
cooper pairs have zero spin projection. The spin dependent BdG equations are decoupled 
into two independent sets of (two component) equations, one for the spin up electron, spin 
down hole quasiparticle {u^{r),Vj^{r)), and the other for {ui{r),v-^{r)). The corresponding 
BdG equations for spin index s(s) =t (|) or s(s) =| (t), read |^| 



Eus{r) = (Heir) - pU{r))u,{r) + / rfr'A,,(s, x)v^{r') 



(4) 
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where U{r) is the exchange potential, p is 1(— 1) for up(down) spins. Ass{s, x) is the matrix 
element of the pair potential, after a transformation from the position coordinates r, r' to 
the center of mass coordinate x = (r + r')/2 and the relative vector s = r — r'. After 
Fourier transformation the pair potential depends on the related wave vector k and x. In 
the weak coupling limit k is fixed on the Fermi surface {\k\ = kp), and only its direction 6 
is variable. After applying the quasi-classical approximation, i.e. 



Us{r) 


] 




Us{r) 


\ 


y Mr) 


) 




^ Vs{r) 


J 



(5) 



so that the fast oscillating part, of the wave function is divided out, the BdG equations are 
reduced to the Andreev equations 

ETLs{r) = -ih^/mk ■ \7Us{r) + Ass{6, r)vs{r) 

(6) 

Evs{r) = ih^/mk ■ syvsir) + A^^{9, r)us{r) 

where the quantities Us{r) and v-s{r) are electron- like and hole-like quasiparticles with spin 
index s, and s respectively. 

We consider the ferromagnet / insulator / superconductor junction shown in Fig. |l|. The 
geometry of the problem has the following limitations. The particles move in the xy-plane 
and the boundary between the ferromagnet (x < 0) and superconductor (x > 0) is the yz- 
plane at x = 0. The insulator is modeled by a delta function, located at x = 0, of the form 
V6{x). The temperature is fixed to K. We take both the pair potential and the exchange 
energy as a step function i.e. Ass{0,r) = Q{x)Ass{0), Ur) = 0(— x)t/. For the geometry 
shown in Fig. [1|, Eqs. ^ take the form 

Eus{x) = -ih^/mkFx£us{x) + Ass{0)vs{x) 
Evsix) = ih^/mkF,£v^{x) + A*-^{9)us{x) 

When a beam of electrons is incident from the ferromagnet to the insulator, with an angle 
9, the general solution of Eqs. (^, is the two component wave function = (miji], t>|[|]) 
which for x < is written as 
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X COS Oa 



+ b 



V) 



(8) 



where a|[j^], are the amphtudes for Andreev and normal reflection for spin up (down) 
quasiparticles, and ^^[j^] = ^^{Ep ± U) is the wave vector of quasiparticles in the ferro- 
magnet for up (down) spin. The wave vector of the electron-hke, hole-hke quasiparticles is 



approximated by kg — y "^'' ■ Since the translational symmetry holds in the y-axis direc- 
tion, the momenta parallel to the interface is conserved, i.e. sin 6* = sin6'A = A;sSin6's. 
Note that 9 is different than 9a since the retroreflection of the Andreev reflection is bro- 
ken. Using the matching conditions of the wave function at x = 0, ^/(O) — ^'//(O) and 
^//(O) ~ ^/(O) — (2my/^^)^'/(0), the Andreev and normal reflection amplitudes am],6m] 
for the spin up (down) quasiparticles are obtained 

4n+Ai 



(-1 - Ai - iz^ii]){-l -X2 + iz-iii]) + (1 - Ai - iz^ii]){-l + A2 - iz^ii])n+n^(f)^(f)*+ ' 



(9) 



(-1 - A2 + izm]){l - Ai + iz^[i]) + (-1 + A2 - izm]){-l - Ai + izt^])n+n_0_0; 
(-1 - Ai - izm]){-l -X2 + iz^ii]) + (1 - Ai - izm]){-l + A2 - i2;m])n+n_0_(/):^ ' 

(10) 



where zq = ^T[il 
given by 



-2a_ Ai = A2 = S2^^pl. The BCS coherence factors are 

cos Os ' COS O3 Ks ' COS Os K3 



u 



± 



= [1 + 4W^\Kd&)?/m. 



(11) 



[1 - .jE^-\^^m'/m. 



(12) 



and n± = v±/u±. The internal phase coming from the energy gap is given by 0± = 
[A±(^)/|A±(6')|], where A+(6') = A(^) (A_(^) = A(7r - 6*)), is the pair potential experi- 
enced by the transmitted electron-like (hole-like) quasiparticle. A(^) = A-fj^(^) = Aj^-f(^), 
since the cooper pairs have zero spin projection i.e. d \\ z. 
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When 6 > sin~^(^) = Od total reflection occurs and the spin and charge current van- 
ishes. When 9ci > 9 > sin^^(^) = 9c2, although the transmitted quasiparticles in the 
superconductor do propagate, the Andreev reflected quasiparticles do not propagate. This 
process is called virtual Andreev reflection (VAR process). In this case the spin and charge 



current do not vanish since a finite amplitude of the Andreev reflection still exists. |T5| 
According to the BTK formula the conductance for the charge current of the junction, 
9), for up(down) spin quasiparticles, is expressed in terms of the probability ampli- 
tudes a^i],b^[i] as 



a,^^^^{E,9)=Re 



1 + T^I«TU]P - I^TUll 



Ai 

The tunneling conductance, normalized by that in the normal state is given by 



(13) 



a,{E) = a,JE) + a,XE), (14) 



^'?T[ii(^) = / d9 COS 9ag {E, 9) P^i^q^i], (15) 

J — 7r/2 



where 



= d9 cos 9[aN,{9)P^q^ + a^^(0)P^gJ, (16) 

J — 7r/2 



^^™^^^ = (i + Ao^'+4,' ^^^^ 

where P^yi] = {Ep ± U)/2Ef is the polarization for up(down) spin. In the zq = limit the 
interface is regarded as a weak link, showing metallic behavior while for large zq values the 
interface becomes insulating. 

III. POSSIBLE SPIN TRIPLET PAIRING STATES 

For the spin triplet pairing state the Cooper pairs have spin 1 degree of freedom. The 
gap function is a 2 x 2 symmetric matrix which in the spin space can be written as 
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A{k) =iay{d{k) ■ &), (18) 

where <t denotes the Pauh matrices and d{k) is a vectorial function which is odd in k. The 
d vector defines the axis along which the Cooper pairs have zero spin projection. In the 
following we will take d \\ z. In that case A-fi = = 0, while A^j = A^^ = A(^^). The 
energy spectrum of the quasiparticles consist of two branches which are identical for unitary 
pairing states, i.e. A\k)A.{k) is proportional to the unit matrix, and distinct for non 
unitary states. The non-unitary states have been ruled out for Sr2Ru04 by the vary small 
residual value of the specific heat at zero temperature [0]. In this paper we will examine 
only the case of unitary pairing states. As an example we consider the state 



A{9) = Ao 



(19) 



a) For the isotropic p-wave pairing state A|j(^^) = Aj|(6') = Aoexp{i{9 — /?)), and 
^Tt(^) ~ ^ii(^) = 0, /5 denotes the angle between the normal to the interface and the 
X-axis of the crystal. This is an opposite spin pairing state, with a gap of constant modulus 
for both spin parts on the Fermi surface. In the following we will consider the cases where 
the matrix element A|j(^^) (expressed as A (6')), of Eq. ^ has the following 6 dependences. 

b) In case of a p-wave superconductor, proposed by K. Miyake, and O. Narikiyo 



A(^) = — [sin(A;^a) +isin(A;j^a)], (20) 

with k^a = Rtt cos{6 — (3), and kya = RTrsm(9 — (3), sm = v^siu;^ = 1.125, and R = 0.9. 
This state does not has nodes. 

We consider also three pairing symmetries for Sr2Ru04 with line nodes. 

c) In the first 2D /-wave state Big x 

A{e) = Aocos2{e - f3)[cos{e - f3) +tsm{e - (3)]. (21) 

This state has nodes at the same points as in the dx^_y2-wa.ve case. 

d) For the second 2D /-wave state i?2g x 
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A{9) = Aosm2{9 - P)[cos{9 - P) + ism{9 - P)]. (22) 



This state has nodes at 0, 11/2, vr, 37r/2, and has also been studied by Graf and Balatsky 
e) In case of a nodal p-wave superconductor 

A(^^) = — -[sm{kxa) +ism{kya)], (23) 
Sm 

with k^a = TT cos{9 — f3), and kya = 7rsin(^ — f3). We use here the same normalization 
proposed by Dahm et al sj\/ = sin ^ = 1.125, where the Fermi wave vector is chosen 



as kpa = IT, in order to have a node in A(^^). This state has nodes as in the x state. 
The corresponding nodeless form was initially proposed by K. Miyake, and O. Narikiyo ||2l[ 
and is considered ClS db separate case. 

IV. TUNNELING CONDUCTANCE CHARACTERISTICS 

In Figs. - ^ we plot the tunneling conductance (Jq{E) for different values of the exchange 
interaction x = U/Ep (a) zq = 0, [3 = 0, (b) zq = 2.5, (3 = 0, (c) zq = 2.5, (3 = tc/A. The 
pairing symmetry of the superconductor is Big x Eu in Fig. 0, x Eu in Fig. |^, nodal 
p-wave, in Fig. ^ p-wave, proposed by K. Miyake, and O. Narikiyo ^T|] in Fig. ^, and 



isotropic p-wave in Fig. ^ When the ferromagnet is a normal metal i.e. x = the results 
of ||2^ are reproduced. For zq = 0, the subgap conductance is suppressed, with the increase 
of X, as in the case of a c?a.2_y2-wave superconductor. |jT5[ 



In the case of normal metal / insulator / triplet superconductor junction the peaks inside 
the gap, are connected to bound states, which are formed due to the sign change that the 
transmitted quasiparticles feel, for fixed (3 at discrete values of 9. The conductance peaks 
occurs at these energies where an increased number of bound states is formed. [^. For 
unitary pairing states the spins of the incident electron and the Andreev reflected hole are 
opposite and since the spin up and spin down quasiparticles have equal wave vectors, no 
spin effects are involved in the Andreev reflection. This is not true, when the normal metal 
is replaced by a ferromagnet. In that case the spin- up and spin-down wave vectors are 
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not equal and the spin affects tlie Andreev reflection. The Andreev reflected hole decays 
in the ferromagnet and the interference with the reflected electron is weak. Moreover the 
transmitted quasiparticles experience weakly the sign change of the pair potential which 
is the reason for the formation of the conductance peak. Due to this the conductance 
peaks are suppressed. This is seen in Fig. |^ b(c), zq = 2.5 and (3 = 0(7r/4), for the 
Big X Eu case. Quantitatively the suppression of the conductance peaks in the ferromagnet 
/ insulator / triplet superconductor junction can be seen if we calculate the magnitude of 
the Andreev reflected amplitude as a function of the exchange field when a bound state is 
formed. Then the amplitude would decay to zero with the increase of the exchange field. 
This calculation has been done in the case of a ferromagnet / insulator / time reversal 
symmetry broken superconductor junction where simple arguments have been derived to 
connect the suppression of the Andreev reflection as x increases with the reduction of the 



conductance peaks. [27 



Also in the metallic limit (x = 0) for all the pairing states we see the presence of a large 
residual density of states within the energy gap as a signature of unconventional pairing 
symmetry with higher than two angular momentum. This is modified by the presence of 
the ferromagnet, where the increase of the exchange field suppresses the density of states 
within the gap. Also in the metallic limit the conductance increases linearly with E, which 
is consistent with the presence on line nodes in the pairing potential. This linear form of the 
spectra remains unchanged when increasing x. Generally the evolution of the conductance 
spectra with the increase of x, for the three pairing symmetries with line nodes depends 
strongly on the position of the nodes in the pairing potential. In the (ij.2„j^2-wave case 
the peak at E = and (3 = 7r/4, due to the sign change of the pairing potential for the 
transmitted quasiparticles, is largely reduced with the increase of x [|I5|. In the case of 
Big X Eu-wave state for a; = 0, the pairing potential is more complicated and the sign 
change occurs at discreet values of 6, for fixed (3. As a result bound states are formed within 
the gap, and also the position of the conductance peaks depends on the orientation angle 
13, as seen in Fig. Also the spectra for angle (3 in the Big x Eu case is identical to the 
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spectra of i?2g x Eu in Fig. ^ for angle n/A — f3, since the nodes for the two symmetries 
differ by tt/A. The evolution with x of the conductance spectra for zq = 0, is different in the 
Big X Eu, X Eu cases. It develops a dip at = in the Big x Eu-, while it has a peak in 
the B2g X Eu case, at = 0. The nodal p-wave case in Fig. ^ has the same nodal structure 
as the B2g x Eu case and we see that the spectra for these two candidates are similar. 

For nodeless pairing states a subdap or a full gap opens in the tunneling spectra. This is 
seen in Fig. H, for the pairing state proposed by K. Miyake, and O. Narikiyo (MN) The 
spectra is similar to the nodal j9-wave case, except that in the (MN) case a subgap opens in 
the tunneling spectra for certain junction orientation. In this region for = the tunneling 
conductance is equal to 2, when x = and has a constant value for x > 0. The tiny subgap 
is an indication of nodeless pairing state. For the isotropic p-wave case the tunneling spectra 
changes with zq, as can be seen in Fig. but not with the boundary orientation (3. The 
spectra is nodeless, and for zq = 0, the conductance is <Jq{E) = 2, within the energy gap, 
for X = 0. Similar results have been obtained in Ref. |T^, where the tunneling conductance 
of a ferromagnet / triplet superconductor interface is calculated, for both unitary and non- 
unitary pairing state, having Eu symmetry. Generally in all pairing states the reduction 
of the subgap conductance with the exchange field is symmetric since the density of states 
modulation within the subgap is not induced by spin dependent effects. 

V. MAGNETIC FIELD EFFECTS 

In this section we describe the effect of the external magnetic field H in the spectra for 
different values of the exchange field x. We will see that since the effect of the magnetic 
field depends on the spin, the evolution of the tunneling spectra with x is asymmetric. The 
tunneling conductance is given by 

ag{E) = a,^ {E - ^ibH) + (Tg, {E + ^ibH). (24) 

In Fig. |^,b,c the tunneling conductance <Jq{E) is plotted for fixed magnetic field ^bH/'^o = 
1, and barrier strength zq = 2.5, for different values of the exchange interaction x. The 
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pairing symmetry of the superconductor is Big x Eu, i?2g x Eu, nodal p-wave, respectively. 
The same information is plotted in Fig. ^ a,b, for the nodeless pairing states, isotropic p- 
wave, andp-wave, proposed by K. Miyake, and O. Narikiyo, |21[| respectively. The orientation 
of the superconductor is chosen as /3 = 0. In the absence of the exchange interaction 
(x = 0) the magnetic field splits symmetrically the tunneling spectrum. The amplitude of 
the spiting depends linearly on the magnetic field H. The main effect of the polarization is 
the imbalance in the peak heights for E positive and negative. The ratio of the peaks for 
positive and negative energy is proportional to the exchange field of the material. For E < 
the pattern is suppressed linearly with the increase of the x, while for E > the tunneling 
conductance spectra initially increases with x and then decreases. Also the conclusions of 
the previous section for the nodal (nodeless) form of the tunneling spectra are still valid in 
the presence of a magnetic field. 



VI. CONCLUSIONS 

We calculated the tunneling conductance in ferromagnet / insulator / triplet supercon- 
ductor junctions, using the BTK formalism. We assumed pairing potentials with nodes such 
as the nodal p-wave and two 2D /-wave states with line nodes, breaking the time rever- 
sal symmetry. Also we examined two nodeless pairing states, the p-wave proposed by K. 
Miyake, and O. Narikiyo and the isotropic p-wave. The linear variation of the conduc- 



tance with E is an indication of line nodes and is not influenced much when the exchange 
interaction increases. On the other hand the large residual density of states within the gap 
is reduced with the increase of x, and the peaks due to the formation of bound states are 
removed due to the suppression of the Andreev reflection. The evolution of the spectra with 
X depends on the position of the nodes and is different in the three pairing states with line 
nodes. In the case of nodeless pairing states the tunneling conductance develops a subgap 
or a full gap where (Jq{E) has a constant value within. The exchange interaction suppresses 
the conductance within the gap, and can be considered as a measure of the polarization of 
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the material. The magnetic field splits linearly the tunneling spectra and in the half metallic 
ferromagnetic limit x — 1, eliminates the negative branch of the spectrum. These features 
can be used to distinguish between the candidate pairing symmetry states of Sr2Ru04. 
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FIG. 1. The geometry of the ferromagnet / insulator / triplet superconductor interface. The 

pairing state is unitary with zero spin projection. The vertical line along the y-axis represents the 
insulator. The arrows illustrate the transmition and reflection processes at the interface. 9 is the 
angle of the incident electron and the normal, 6a is the angle of the reflected hole and the normal, 
and 9s is the angle of the transmitted quasiparticle and the normal. Note that 9 is not equal to 9a 
since the retroreflection of the Andreev process is lost. 
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E/Ao 

FIG. 2. Normalized tunneling conductance crq{E) as a function of E/Aq for x = (solid 
line), X = 0.4 (dotted line), x = 0.8 (dashed line), and x = 0.999 (long dashed line), for different 
orientations (a) Z=0, (3 = 0, {h)Z = 2.5, (3 = 0, (c) Z = 2.5, f3 = ir/A. The pairing symmetry of 
the superconductor is Big x E^- 
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FIG. 6. Normalized tunneling conductance CTq{E) as a function of E/IS.Q for x = {) (solid 
line), X = 0.4 (dotted line), x = 0.8 (dashed line), and x = 0.999 (long dashed line), for different 
orientations (a) = 0, /? = 0, (b)2;o = 2.5, /3 = 0. The pairing symmetry of the superconductor is 
isotropic p-wave. 
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FIG. 7. Normalized tunneling conductance (Tq{E) as a function of E/Aq for x = (solid line), 
X = 0.2 (dotted line), x = 0.4 (dashed line), and x = 0.999 (long dashed line), for zq = 2.5, f3 = 0.0, 
for different nodal pairing states (a) Big x Eu, (b) B2g x Eu, (c) nodal p-wave. 
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FIG. 8. Normalized tunneling conductance CFq{E) as a function of E/Aq for x = (solid line), 
X = 0.2 (dotted line), x = 0.4 (dashed line), and x = 0.999 (long dashed line), for zq = 2.5, (3 = 0.0, 
for different nodeless pairing states (a) isotropic p-wave, (b) nodeless p-wave proposed by M.N. 
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